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Abstract. Motivated by recent experimental measurements on pressure-driven phase 
transitions in Mott-insulators as well as the new iron pnictide superconductors, we 
show that first principles Car-Parrinello molecular dynamics calculations are a powerful 
method to describe the microscopic origin of such transitions. We present results for 
(i) the pressure-induced insulator to metal phase transition in the prototypical Mott 
insulator TiOCl as well as (ii) the pressure-induced structural and magnetic phase 
transitions in the family of correlated metals AFe2As2 (A=Ca,Sr,Ba). Comparison of 
our predictions with existing experimental results yields very good agreement. 
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1. Introduction 

Understanding the microscopic origin of phase transitions in correlated systems is a 
challenging task both for theory and experiment due to the large number of degrees of 
freedom (electronic, lattice and orbital) involved in the transitions. In a Mott insulator 
a transition from an insulating to a metallic or even a superconducting state can be 
induced by doping the system with carriers as is the case of the high-T c cuprates pQ, 
or by applying pressure as in the K-charge transfer salts [2J [3] . Pressure-driven Mott- 
insulator to metal transitions have also been discussed in a few transition metal oxides 
like MnO [I], V2O3 [5], and lately in TiOCl [El [7j. Very recently, high temperature 
superconductivity has been achieved in a new class of correlated metals, the iron pnictide 
compounds, both by doping as well as by setting the systems under pressure [El [9]. The 
origin of this phase transition is still under strong debate. 

External pressure- induced phase transitions are especially attractive for microscopic 
modelling since the chemistry of the system remains untouched, i.e. possible disorder 
as in the doping case is almost inexistent, and one can concentrate on the lattice 
changes induced by pressure. A reliable microscopic description of such transitions 
requires a simultaneous treatment of the electronic and lattice degrees of freedom as is 
realized in the first principles Car-Parrinello molecular dynamics method [10J. Instead of 
considering the motion of the nuclei and the solution of the Kohn-Sham equation for the 
electrons at fixed atomic positions as separated problems, the Car-Parrinello method 
explicitly views the electronic states as dynamical variables and considers these as a 
unified problem. By writing a fictitious Lagrangian for the system which leads to coupled 
equations of motion for both ions and electrons, the lattice dynamics and electronic 
properties can be treated on the same footing. Furthermore, the Parrinello Rahman 
approach to molecular dynamics simulations at constant pressure [11] treats the lattice 
parameters as additional degrees of freedom, allowing the first principle observation of 
structural changes as a function of pressure. The great advantage of this procedure is the 
ability to deal with crystals with anisotropic compressibilities along different directions 
as is the case of TiOCl and the new iron pnictide superconductors. In the present 
work we will show by reviewing recent calculations on pressure-driven transitions in 
the multiorbital Mott insulator TiOCl as well as the family of iron pnictides v4Fe2As2 
(A=Ca,Sr,Ba) that the Car-Parrinello method is often superior to the more common 
Born-Oppenheimer methods for describing the features of the phase transitions. 

2. Overview on TiOCl and AFe 2 As 2 

TiOCl and AFe2As2 (A=Ca,Sr,Ba) are both layered compounds containing 3d transition 
metal ions as shown in Fig.[H Correlation among the open d shell electrons is responsible 
for a rather unconventional behavior in these materials P, El [121 E2 El EE EE CEO EE 
[111 [201 [HI 1221 [231 [Ml [251 [261 [23 In TiOCl, 

titanium (Ti) and oxygen (O) form bilayers stacked along the c direction and separated 



Pressure-driven phase transitions in correlated systems 3 




Figure 1. Lattice structures of (a) TiOCl and (b) AFe2As2 (y4=Ca,Sr,Ba) for the 
example of BaFe2As2- 

by chlorine (CI) layers, while in v4Fe2As2 (^4=Ca,Sr,Ba), layers of Fe-As tetrahedra are 
separated by alkaline earth metal ions. 

2.1. TiOCl 

TiOCl is a quasi-one- dimensional Mott insulator at room temperature where Ti chains 
along the b direction are weakly coupled [121 H3]- Upon lowering the temperature, 
TiOCl undergoes two consecutive phase transitions from a Mott insulating state to 
a spin-Peierls phase with dimerized Ti chains through an intermediate structural 
incommensurate phase with dominant out-of-chain incommensurability even under 
zero external magnetic field [HI [15J, [HI [17]. Since the discovery of this unusual 
intermediate phase, many efforts have been devoted to understanding its microscopic 
origin jlUJ, HU H2J H3j. Recently, a two-dimensional frustrated spin-Peierls model was 
proposed, where the model parameters were calculated from ab initio density functional 
theory (DFT) [ID] . Besides the dominating interactions along chains, it was pointed out 
that the interchain magnetic frustrations may play an important role in the formation of 
the unconventional incommensurate phase in addition to the interchain elastic couplings. 

Due to the 3c? 1 configuration of Ti 3+ and quenching of the orbital degrees of 
freedom [161 E]> TiOCl was thought to be a model system for the one band Hubbard 
model. In order to investigate a possible Mott insulator-to-metal transition in this 
system in the framework of the one-band Hubbard model various experimental groups 
have recently investigated the effect of pressure [6l [7j EUJ EH [22] and electron 
doping [IH1 [19] on TiOCl. Interestingly, pressure driven experiments lead to different 
conclusions. Kuntscher et al. [HI [22] reported from optical measurements and X- 
ray diffraction analysis the observation of an insulator-to-metal transition at 16 GPa 
accompanied by a structural symmetry change. Forthaus et al. [7J, on the other 
hand, neither observed any indication of a metallic state up to pressures of 24 GPa 
from transport measurements nor did they identify from X-ray diffraction analysis 
changes on the crystal symmetry at 16 GPa. Recently, Blanco-Canosa et al. [21] 
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found a dimerized insulating state at high pressure by X-ray diffraction analysis 
and magnetization measurements combined with ab initio DFT results. The DFT 
calculations in Ref. [21] consisted in performing a structural optimization within the 
generalized gradient approximation (GGA) and an analyis of the electronic properties 
within GGA+U. It should be stressed that GGA gives a metallic state for TiOCl while 
GGA+U describes TiOCl as an insulator. In the following, we will argue that inclusion 
of strong correlation is crucial for obtaining the correct optimized lattice structure. 
Therefore, the theoretical understanding of the dimerization under pressure |21j may be 
questionable due to the alternate use of different functionals in one study. 

In a recent work [20], we resolved the previous inconsistencies by performing ab 
initio Car-Parrinello molecular dynamics [TU], [TTJ within the projector-augmented wave 
basis (CP-PAW) [45] on TiOCl under pressure. The GGA+U functional was considered 
throughout the whole study. Within GGA+U correlation effects are included at the 
mean-field level. The simulations show that two phase transitions occur; a first one 
from the Mott insulating state to an intermediate dimerized metallic state accompanied 
by a symmetry lowering from orthorhombic Pmmn to monoclinic P2i/m, and a second 
transition to a normal metal with recovered Pmmn symmetry. The first phase transition 
is consistent with the observations of Kuntscher et al. [H [22] and the structural phase 
transition is in agreement with the experimental result of Ref. [21] . The second phase 
transition has not yet been reported experimentally. Moreover, the pressure effects 
involve the three t 2g orbitals in TiOCl. 

2.2. AFe 2 As 2 

The materials AFe 2 As 2 (A=Ca,Sr,Ba) are paramagnetic metals at room temperature. 
With decreasing temperature, a magnetic phase transition to a stripe-type 
antiferromagnetic (AF) ordered state and a structural phase transition from tetragonal 
to orthorhombic symmetry take place simultaneously [231 1211 [25], [26] . Even more exciting 
is the fact that these systems are superconductors both under pressure j9[ [3U [35[ [361 
1371 [38] and electron or hole doping [281 [211 [3Ql [31], which has attracted considerable 
and persisting interest. 

Since doping usually induces disorder, it is more practical to investigate the origin 
of the phase transition and the phenomenon of superconductivity by applying external 
pressure. However, while the pressure- driven magnetic and structural phase transitions 
in CaFe 2 As 2 have been thoroughly investigated [321 [33] , those in SrFe 2 As 2 and BaFe 2 As 2 
have not yet been fully studied due to the high Sn content that the grown crystals 
have [31]. Therefore, a theoretical prediction of the phase transitions in SrFe 2 As 2 
and BaFe 2 As 2 is desirable since it is free of the Sn-problem. Unfortunately, it is 
well-known that the magnetization and the volume of the iron pnictide compounds 
are overestimated by DFT calculations within GGA when structural optimizations are 
performed [4"6 l I4"T [ l4~8] . We will show here that due to this overestimation we only have to 
simulate higher pressures than the experimental ones in order to reach the experimental 
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volume conditions while structural changes at critical pressures and the main features 
of the electron properties within different phases can still be captured properly. 

Optimization of cell parameters and atomic positions under pressure within the 
framework of density functional theory (DFT) has been done by Yildirim [49] and Xie et 
al. [50] ■ By considering Vanderbilt-type ultrasoft pseudopotentials, a smooth structural 
transition for CaFe 2 As 2 was observed under pressure without detection of a sudden 
sizable increase of the cell parameters a and b or a strong decrease of the cell parameter 
c which is inconsistent with experimental results [32] . Xie et al. [50] optimized within the 
full potential linearized augmented plane wave method (FPLAPW) the orthorhombic 
lattice structure for BaFe2As2 under pressure by relaxing the internal parameter za s 
and the c/a ratio while keeping the b/a ratio fixed. This procedure doesn't allow for 
the detection of the structural and magnetic phase transitions. 

We investigated the phase transitions of ^4Fe2As2 (A=Ca,Sr,Ba) under pressure 
in the lower temperature region with the CP-PAW method within the spin-polarized 
GGA functional. The first order phase transition from an orthorhombic phase 
to a collapsed tetragonal phase in CaFe 2 As2 is confirmed with relative changes of 
lattice parameters, bond lengths and angles agreeing very well with the experimental 
observations, indicating the validity of the application of the CP-PAW method to this 
system. At the critical pressure, an abrupt disappearance of magnetization is also found 
as observed experimentally. In particular, our calculations can account for the sudden 
expansion in the ab plane observed experimentally [32] at elevated pressure crossing the 
critical pressure as we discussed in ref. [47J, which was not obtained in other calculations. 
Applying this method to study SrFe 2 As 2 and BaFe 2 As 2 , where less is known about the 
features of the phase transition, a simultaneous structural (orthorhombic to tetragonal) 
and magnetic phase transition is predicted at high pressure. However, we observe a weak 
first-order phase transition for SrFe2As2 and a continuous phase transition for BaFe2As2 
in contrast to the strongly first order phase transition in CaFe 2 As 2 . Finally, it should 
be pointed out that multi-orbital physics is important in AFe2As2 (A=Ca,Sr,Ba) since 
all five d orbitals cross the Fermi level. 

3. Details of our calculations 

In our studies of both systems, full, unbiased relaxations of all lattice and electronic 
degrees of freedom have been performed at each pressure value. In TiOCl, local 
coordinates are chosen as x = b, y = c, z = a while in AFe 2 As 2 (v4=Ca,Sr,Ba), x — a, 
y = b, z = c. 

All our calculations were performed with time steps of 0.12 fs at zero temperature. 
The system size was 4 x 4 x 4 k points for AFe2As2 (A=Ca,Sr,Ba) and 6x6x6 for 
TiOCl on doubled unit cells corresponding to different AF configurations [201 HZ]- We 
used high energy cutoffs of 612 eV and 2448 eV for the wave functions and charge density 
expansion, respectively. The total energies were converged to less than 0.01 meV/atom 
and the cell parameters to less than 0.0005 A. In TiOCl, a value ofU = 1.65 eV was used 
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Figure 2. Comparison of the density of states for TiOCl between CP-PAW (a) and 
FPLAPW (b) methods within GGA. 



for the GGA+U calculations with which correct spin exchange along b is reproduced [40J. 
In AFe2As2 (A=Ca,Sr,Ba), the 3s3p3d (4s4p4d/5s5p5d) states in Ca (Sr/Ba) are treated 
as valence states and 3d4s4p in Fe and As, while in TiOCl, the 3s3p3d4s in Ti, 2s2p 
in O and 3s3p in CI. We have checked our calculations within GGA based on these 
choices of valence states with the FPLAPW method as implemented in the WIEN2k 
code [51]. Very good agreement was found between these two methods, indicating the 
correctness of our constructed PAW basis for both systems. In Fig. [2} we show one of the 
comparisons for TiOCl between CP-PAW and FPLAPW within GGA. It is found that 
the total density of states (DOS) and partial DOS (pDOS) calculated from FPLAPW 
can be reproduced by CP-PAW. Not only the peak position but also the shape of the 
DOS are almost the same throughout the whole energy range. More spikes in the DOS 
from CP-PAW come from less k points used in the calculation. 

Table 1. Comparison of equilibrium lattice parameters for TiOCl optimized by 
different functionals. B3LYP data were obtained from Ref. [33] 





a (A) 


b(A) 


c(A) 


exp. 


3.79 


3.38 


8.03 


CRYSTAL B3LYP 


3.81 


3.49 


8.69 


GGA 


3.89 


3.30 


8.06 


GGA+U(1.65,FM) 


3.80 


3.49 


8.12 


GGA+U(1.65,AF) 


3.81 


3.42 


7.51 



Now let us investigate the effect of strong correlations on the lattice optimization [52] 
of TiOCl. In Table. [H we present the optimized lattice parameters within different 
functionals compared to the experimental equilibrium state in TiOCl. For the lattice 
parameter c, we find that it is rather hard to determine precisely due to the fact that 
there only exist weak van der Waals forces between layers along c. Thus, by varying 
the lattice parameter c, the total energy remains almost unchanged. For the lattice 
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Table 2. Comparison of equilibrium atomic distances and angles for TiOCl optimized 
by different functionals. B3LYP data obtain from Ref. [33] 

Ti-0| a (A) Ti-Cl (A) Ti-0-Ti| a 



exp. 


1.96 


2.40 


150° 


CRYSTAL B3LYP 


1.98 


2.45 


148.6° 


GGA 


1.98 


2.48 


156.4° 


GGA+U(1.65,FM) 


1.97 


2.44 


149.1° 


GGA+U(1.65,AF) 


1.97 


2.41 


151.4° 



parameters a and b, B3LYP and GGA+U with both FM and AF configurations can 
provide a better comparison to experiment than GGA does, indicating the importance 
of strong correlations which are absent in GGA. Futhermore, it is found that the results 
from GGA+U with an underlying AF ordering is even better than that from GGA+U 
with FM ordering, which indicates that the correct setting of the spin configuration 
improves the lattice optimization. In Table. [2} we show the optimized atomic distances 
and angles by different functionals compared to the experimental results. Again, the 
same trend is found as in Table. [U i.e., with the inclusion of correlations as implemented 
in GGA+U and in B3LYP, the optimized structures are closer to the experimental one. 

In contrast to TiOCl, the iron pnictide compounds are metallic even in the AF 
state and are considered as moderately correlated systems. Therefore, the lattice 
optimizations for AFe2As2 (A=Ca,Sr,Ba) have been performed within GGA. 

4. Results 

4-1. TiOCl under pressure 

In the study of TiOCl under pressure, we use a renormalized pressure due to the 
overestimation of the critical value. However, it has a different origin compared to that in 
AFe2As2 (A=Ca,Sr,Ba). As reported from experiments [13], TiOCl is a Mott insulator 
at room temperature and ambient pressure. Only short-range AF spin fluctuations 
exist in TiOCl. However, to reach the insulating state by DFT calculations, one has to 
impose certain long-range spin ordering. It is well-known [53] that the AF insulator is 
more robust than the paramagnetic Mott insulator under bandwidth control; one has 
to apply much higher pressure to induce an insulator-to-metal transition if long-range 
AF ordering is present. 

In Fig. [3] we present the calculated change of atomic distances as a function of 
applied external pressure. The lattice parameters at each pressure are obtained by fully 
relaxing the lattice structure without any symmetry constraints. Although a magnetic 
ordering has to be imposed within the GGA+U functional (AF ordering along b and 
FM configuration along the other directions) the results are completely consistent with 
the experiments (TJ [22] . 

At P = P c (black dashed vertical line in Fig. [3]), the equivalent Ti-Cl distances in 
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Figure 3. Calculated change of atomic distances of TiOCl as a function of applied 
external pressure, normalized to the critical pressure, (a) The distance between Ti and 
CI. (b) The distance between Ti and O along a. 
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Figure 4. Density of states of TiOCl in three different phases, (a) is for the 
low pressure Mott insulating state (P < P c ), (b) for the dimerized metallic phase 
{P c < P < P' c ), (c) for the uniform metal (P > P' c ). 



the lower pressure region split into four different values (see left panel) and the Ti-0 
distance into two (see right panel), indicating a lowering of the symmetry. We found 
that the space group changes from Pmmn to P2\/m which agrees with the experimental 
observation [21]. Most interestingly, both inequivalent Ti-0 distances show jumps at 
this phase transition, indicating that the system expands along a in agreement with 
experiment [54J. Further increasing the external pressure, we find a second jump of the 
atomic distances at P = P' c > P c (blue dotted vertical line in Fig. [3]). At this pressure, 
all the inequivalent Ti-Cl and Ti-0 distances become equivalent again, indicating the 
recovery of the Pmmn symmetry. This phase transition has not yet been observed. 

Fig. Hj) presents the DOS in these three phases. Fig. Hj) (a) shows the DOS for the 
lower pressure case. The DOS looks very similar to that at ambient pressure except for 
the widening of the bands and the reduction of the gap amplitude. Right above the first 
phase transition, the gap is closed as shown in Fig. HJ) (b), which means that the system 
should show metallicity. The obtained insulator-to-metal transition is consistent with 
optical experiments P, [22] but contradicts the electrical resistivity measurements [7]. 
A possible reason is that TiOCl is a strongly anisotropic crystal with a large van der 
Waals gap along c. Thus, if the resistivity measurements have been done along c, the 
system should show insulating behavior, reminiscent of cuprates. Further increasing the 




Figure 5. Comparison of DOS for CaFe2As2 calculated for optimized and 
experimental lattice structures in both AF orthorhombic phase ((a), (b)) and 
paramagnetic 'collapsed' tetragonal phase ((c), (d)) at low temperature. More detailed 
comparison of the pDOS of Fe d states is shown in (e) and (f). 



pressure to P > P' c , it is shown in Fig. H]) (c) that the metallicity remains with further 
widening of the bands and even nearly closure of the gap between the p states of O, CI 
and the d states of Ti. Therefore, this phase transition is exclusively induced by the 
structural change. 



4-2. AFe2As2 underpressure 

In Fig. [5j we show the comparison of densities of states calculated based on optimized 
and experimental lattice structures for CaFe2As2 in both the AF orthorhombic phase 
and the paramagnetic 'collapsed' tetragonal phase at low temperature. Figs. [5] (a) and 
(b) show comparison of DOS in the orthorhombic phase. It is found that the DOS are 
almost the same although the volume and the magnetic moment are both overestimated 
by DFT calculations within the spin-polarized GGA functional. Figs. [5] (c), (d) are 
the comparison of DOS in the 'collapsed' tetragonal phase, and Figs. [5] (e), (f) the 
corresponding comparison of pDOS of Fe d states. It is obvious that both DOS and 
pDOS are almost identical. Peaks coming from Fe d xz /d yz and d x 2_ y 2 are shifted 
away from the Fermi level due to the structural change under pressure which avoids 
the high instability at the Fermi level present in the high-temperature tetragonal 
phase. Considering the fact that the change of the lattice structure under elevated 
pressure observed experimentally can be described well by our CP-PAW calculations, in 
particular the jumps of the volume and the distances of Fe-As at the phase transition 
which are the two key quantities determining the volume collapse and the electronic 
properties, respectively, (as shown in Ref. [S]), we come to the conclusion that our 
DFT structural optimizations are valid in both low temperature phases, namely, the 




Figure 6. Calculated changes of (a) volume and (b) magnetization of AFe2As2 
(^4=Ca,Sr,Ba) as a function of pressure, normalized to the corresponding critical 
pressures. 



AF orthorhombic phase and the 'collapsed' tetragonal phases, except that we have to 
apply higher pressure than the critical value obtained from experiments to eliminate 
the overestimated magnetic moment and volume at ambient pressure and to reach 
the experimental volume conditions. In order to avoid confusion, we renormalized the 
pressure with respect to the critical value throughout this work. 

In Figs. [6] (a) and (b) we present the calculated changes of the volume and 
magnetization of AFe2As2 (v4=Ca,Sr,Ba) as a function of pressure. In the case 
of CaFe2As2, the volume and the magnetization decrease gradually with increasing 
pressure and show a discontinuity at the critical pressure, where the system undergoes 
simultaneously the structural and magnetic phase transitions from an AF orthorhombic 
phase to a volume 'collapsed' paramagnetic phase. Our results are in very good 
agreement with experimental data [32] with a volume collapse of AV th ~ 4.1%, 

In Fig. El we also show the predictions for SrFe2As2 and BaFe2As2. We find that the 
structure and the magnetic phase transitions under pressure still occur simultaneously 
in these two compounds, but the nature of the phases is distinctly different from the 
CaFe2As2 case. This has not yet been confirmed experimentally due to the problem of 
high Sn content [31] in the samples. SrFe2As2 shows smaller magnetization and volume 
jumps at the critical pressure compared to CaFe2As2 while we can hardly detect any 
discontinuity in BaFe2As2, which in this respect resembles more closely jH] the parent 
compound of the iron pnictide superconductors, LaFeAsO. 

In summary, by considering ab initio Car Parrinello molecular dynamics 
calculations, we were able to describe the microscopic details of the pressure-induced 
phase transitions in TiOCl and the family AFe2As2 (A=Ca,Sr,Ba) with very good 
agreement to experimental observations. The CP-PAW proves to be a reliable and 
powerful method to describe structural and magnetic changes at the phase transition. 
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